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We derive a kinetic equation with a non-Markovian collision term which includes a memory 
effect, from Kadanoff-Baym equations in 4 theory within the three-loop level for the two-particle 
irreducible (2PI) effective action. The memory effect is incorporated into the kinetic equation by a 
generalized Kadanoff-Baym ansatz. Based on the kinetic equations with and without the memory 
effect, we investigate an influence of this effect on decay of a single particle excitation with zero 
momentum in 3+1 dimensions and the spatially homogeneous case. Numerical results show that, 
while the time evolution of the zero mode is completely unaffected by the memory effect due to a 
separation of scales in the weak coupling regime, this effect leads first to faster relaxation than the 
case without it and then to slower relaxation as the coupling constant increases. 
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The vacuum of the quantum chromodynamics (QCD) is believed to undergo a phase transition from the hadronic 
• phase to the quark-gluon plasma (QGP) at high temperature T and/or at high quark chemical potential \x, and such 
\ a new state of matter is expected to be produced in heavy-ion collision experiments ongoing at Relativistic Heavy-Ion 
Collider (RHIC) and in the future at Large Hadron Collider (LHC) [1]. Motivated by these experiments as well as the 
inflationary dynamics in the early universe [2-6] , quantum field theories in and out of thermal equilibrium have been 
extensively studied. Many theoretical explanations related to the experiments have been made on the basis of thermal 
and chemical equilibrium. Here, there is an open question about justifying that QGP produced in a transient state 
achieve the equilibration before the phase transition to the hadronic matter. Thermalization of QGP from the initial 
state away from thermal equilibrium [7] has been studied on the basis of the classical Boltzmann equations [8-11]. 
However, it has not been established yet. It might be important to understand the early stages in heavy-ion collisions 
■^J" i from the aspect of nonequilibrium phenomena, compared to a thermally equilibrated state having no information 
' about the initial conditions. 

. The classical Boltzmann equations are often used for studying nonequilibrium properties. Their collision term is 
O ' Markovian one which treats scatterings between particles as independent ones, and no past information is needed to 
calculate it. As a consequence, they cannot be applicable to the strongly correlated or dense plasmas. Undoubtedly, 
generalized kinetic equations are necessary in order to overcome their limitations. 

The quantum Kadanoff-Baym (KB) equations describe the time evolution of the two-time Green's functions which 
involve the information about both the dynamical spectral function and the one-particle distribution function in the 
correlated quantum plasmas [12], and are designed to study time-dependent nonequilibrium phenomena [13-26]. In 
general, it is difficult to solve the KB equations directly due to two-time correlation. It is easier to analyze kinetic 
equations for Wigner distribution function N which can be derived from the KB equations in the equal-time limit, 
by imposing some approximations. When deriving the kinetic equations, one encounters the reconstruction problem: 
How can two-point Green's functions be expressed as a function of N ? While the simplest solution is the Kadanoff- 
Baym (KB) ansatz which is mostly used for deriving Boltzmann-like equations, we have a more general solution, 
the generalized Kadanoff-Baym (GKB) ansatz proposed by Lipavsky et al. [27]. Although the KB ansatz leads to 
inconsistencies in time arguments of a memory effect in the collision term, the GKB ansatz is able to take into account 
it properly, as we will see below. This effect is expected to play an important role in correlated quantum systems 
and/or short-time phenomena. 

The purpose of this paper is to study an influence of the memory effect on the time evolution of the distribution 
function, and to point out that relaxation is strongly affected by the this effect, on the the basis of the kinetic 
equations. Our starting point is a simple question: " If successive collisions among particles in a system of interest 
are not independent ones unlike in the Boltzmann equations, how does relaxation of the system be affected ? " A 
non-Markovian extension of the Boltzmann collision term is derived by using the GKB ansatz, and one has to perform 
a temporal integration over the distribution function having past times in it, as in the original KB equations. However, 
computational efforts are reduced remarkably compared to the KB equations because the distribution function has 
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to be known and time evolved only along the one-time direction. On the other hand, in the KB equations, Green's 
functions have two time arguments and time evolved in the two-time plane. Usually, one expects that, in the weak 
coupling regime the relaxation time can be evaluated correctly by the Boltzmann equations due to a separation of 
scales. However, as we will see, there is the onset of the memory effect in the region where the separation of scales is 
satisfied. 

The present paper is organized as follows. In Sec. II, the KB equations for 4 theory within the 3-loop level for 2PI 
effective action, and the KB and GKB ansatz as approximate solutions of the reconstruction problem arc reviewed. 
Then, with the GKB ansatz, non-Markovian and Markovian collision terms are constructed, and kinetic equations 
with and without the memory effect are derived. Further, those are reduced to linear equations for the distribution 
function in the case of a single particle excitation. Section III deals with numerical calculations. After renormalization 
by the modified MS scheme in a gap equation, the relation between the thermal mass and the coupling constant is 
determined sclf-consistcntly by solving the gap equation numerically. With the solutions of the gap equation, the 
time-dependent damping rate with zero momentum which is a key quantity in the kinetic equations is calculated for 
several values of the coupling constant. Then, the time evolutions of the zero mode are obtained from the kinetic 
equations with/without the memory effect. By changing the coupling constant, an influence of the memory effect on 
the relaxation of the excitation is examined. Some discussions about the numerical results are given at the end of this 
section. Section IV is devoted to summary and outlook. 



II. FROM KADANOFF-BAYM EQUATIONS TO KINETIC EQUATIONS 

KB equations are designed to study noncquilibrium phenomena of a quantum system in terms of one-body Green's 
functions. In this section, we study how kinetic equations are derived from KB equations, by imposing some approx- 
imations. We will see that the memory effect which is neglected in the classical Boltzmann equations is incorporated 
into kinetic equations by the GKB ansatz. 



A. KB equations for 4 theory 

We consider a simple scalar field theory with a Lagrangian density in 3+1 space-time dimensions, 

C = \d^{x)d^{x) ^V(X) + A 4 (x) {21) 

where x = (xq,x), to is the bare mass and A denotes the bare coupling strength. This is an extensive example of a 
relativistic quantum field theory which allows us to challenge theoretical calculations without encountering troubles 
of gauge invariance. 

Within the three-loop order for 2PI effective action, KB equation for G < on Schwinger-Keldysh contour [12,28,29] 
takes a following form in the spatially homogeneous case [13,19,30], 

rxo 

(d 2 XQ +p 2 + m 2 + Z 5 (x ,x )) G < (p,x ,y ) = -i / dz [S > (p, x , z ) - £< (p, x , z )] G < (p,z ,y ) 

Jt 
ryo 

+i dz Y. < (p,xo,z )[G > (p,z ,yo) -G < {p,z ,y )] , (2.2) 
Jt 

where E* 5 (xq, xo) is a tadpole self-energy in the leading order of the coupling constant and £> (p, xo,yo) are sunset 
self-energies in the next-to leading order, and corresponding diagrams are shown in Fig. I. These self-energies are 
expressed in terms of G > as 

A f d 3 k 

T, s (xo,x ) = -J j^G<(k,x ,x ) , (2.3) 
<, \ A 2 f d 3 ki d 3 k 2 d 3 k 3 , N ,,,, Wl 

£> (P,x ,y ) = --J (^(2-)^ (3) (fei + k 2 + k 3 -p) 

xG> (k 1 ,x ,y )G> (k 2 ,x ,y )G> (k 3 ,x ,y ) • (2.4) 
Since G> are full Green's functions, the self-energies are determined self-consistently. We shall discuss this later. 
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FIG. 1. Left and right diagrams represent tadpole and sunset self-energies, respectively. Heavy lines denote dressed propa- 
gators. 

Although one can analyze the KB equations directly, it is easier to consider kinetic equations for the Wigner 
distribution function, which we are concerned with here. For this purpose, it is useful to observe the KB equation 
(2.2) in the equal-time limit x — yo = t, 

(d 2 XQ + p 2 + M 2 (x Q )) G < (p,x a ,y )\ Xo=yo=t = -i / dz a [S>(p, t, z )G < (p, z , t) - £< (p, x , z )G > (p, z , y )] , 

J t 

(2.5) 

where M 2 (x n ) = m 2 + H s (xo, xo)- Considering the similar equation in which the temporal derivative acts on y and 
taking the difference of the two equations, we obtained the following equation [30], 

(d 2 a - d 2 Q ) G < (p,x a ,yv)\ XQ=yo=t = -i / dz Q [S>(p, t, z )G < (p, z Q , t) - £< (p,t, z )G > (p, z ,t) 

J to 

-G > (p,t,z )S < (p,z ,i) +G < (p,t,z )S > (p,z ,i)] . (2.6) 



B. Generalized Kadanoff-Baym ansatz 



In order to derive the kinetic equation satisfied by the Winger distribution function N(p,t) from the KB equation 

(2.6), it is needed to express G > (p, x , y ) as a function of N(p, t). This is the so-called reconstruction problem. The 
simplest solution of this problem is the common Kadanoff-Baym (KB) ansatz [12], given by 



G < (k,x ,y ) = 
G > (k,x ,y ) = 



1 



2w h (t) 
1 

2uJ k (t) 



e -^ k (t) s «N{k, t) + e iw feW S0 (1 + N(-k, t)) 
e -iu k (t)s (1 + N ^ ^ + e ™ k W>oN(-k, t) 



(2.7) 



with so =xq — yo,t= (xq + yo)/2, and Wfc(t) = \J k 2 + M 2 (t). Here, we have adopted the quasiparticle approximation 
and decomposed components of the distribution function with the positive and negative energy. However, this ansatz 
leads to inconsistencies in the time argument of the memory effect because the distribution function depends only on 
the central time t. 

As a more general ansatz, a generalized Kadanoff-Baym (GKB) ansatz has been proposed by Lipavsky et al. [27], 
which is given by 



G < (k,x , j/o ) 



g(fo) 
2w fc (2/o) 

, *(-*(>) 



G < (k,x ,y ) = 



2uj k (x ) 
0(so) 
2^fc(yo) 

2w fc (x ) 



e -*"*(i/o)«o jy^ yo ) + e ^k(vo)so ^ + N(-k, yo)) 

e -<"fc(«o)aoj V ( fe} XQ ) + e <a, fc (so)«o (1 + N(-k, X )) 



-iwu(yo)so 



(1 + N(k, yo)) + e^k(v^ N (-k, y ) 



e -i Wfe (x ) S0 ( X + N ( kj XQ ^ + e *-fc(*o)«o jv(_fe, x ) 



(2.8) 



within the quasiparticle approximation. Note that the GKB ansatz differs from the KB ansatz (2.7) in the time 
argument of distribution functions. This is why the GKB ansatz is able to takes into account the memory effect 
properly. Therefore, we will use the GKB ansatz below. 
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C. Kinetic equation with/without the memory effect 



In the GKB ansatz (2.8), the KB equation in the equal-time limit (2.6) for the positive energy component N(p,t) 
reduces to 



d_ 

dt 



N(p,t) 



6 



d 3 kx d 3 k 2 
(2tt) 3 (2tt) 3 



1 



Vk^p-k.-k^p 



N -kN-k 2 N k 1+ k 2 -p N P 



N_ ki N_ k2 N ki+k2 _ p N p 



: {cos(u fel +uj k2 +uj p - kl -k 2 -u! p ){t-ZQ) 

+ COs(w fcl + UJ k . 2 +UJp-k 1 -k 2 +u) p )(t - z ) 

+3cos(o; fel +uj k2 -^ v - kl - k2 - uj p )(t - z Q ) [N ki N k2 N ki+k2 _ p N p - N ki N ka N ki+k2 _ p N p ] 
+3cos(w fcl - Lu k2 - u v - kl - k2 - u p )(t - z ) [N kl N_ k2 N ki+k2 _ p N p - Wj bi JV_ fe iVj bi+fca _ p JVp] } 



(2-9) 



where ui k = u) k {zo), N k = N(k,z ) and N k = 1 + N(k 7 z a ). The right-hand side of (2.9) is the non-Markovian 
extension of the standard Boltzmann collision term. The memory effect is described by the integration over past 
times which the distribution functions have. This effect is expected to be evidently important if the relaxation time 
of the system of interest is shorter than the memory time itself, and its importance has been illustrated in various 
non-relativistic systems [31,32]. For the KB ansatz, the time argument of distribution functions and energies in the 
integral would have been (zo +t)/2 instead of zq, which is inconsistent with the time argument of the memory effect 
as we mentioned in the introduction. 

If neglecting the memory effect, that is, N(k, z ) — ► N(k, t) and uj k (zo) — ► u> k (t) in the collision term, the integration 
over zq can be performed analytically, and we obtain the finite-time generalization of the Boltzmann equation with 
the Markovian collision term, 



dt 



6ui z 
sin(w fcl + ^fe 2 



d 3 h d 3 k 2 



(27r) 3 (27r) 3 8uJ kl UJ k2 U)p- kl - k 
Wp_fe 1 _fe 2 - UJp)(t - t ) 



+ 



+3 
+3 



Wfe! + ^k 2 + Up- kl - k2 - UJ p 

sin(o; fcl + uj k2 + uj p „ kl ^ k2 + <jj p )(t - t ) 

Wfej + UJ k2 + UJ p - kl - k2 + OJp 

sin(^fci + u k2 - LO p - kl - k2 - LOp){t - t ) 

Wfej + LO k2 — Ul p -k 1 -k 2 — UJ p 

sin(w fcl - ^fe 2 - w P -fc 1 -fc 2 - u P )(t - t ) 



N h N k N n u u N p ~ N k Nu N n u u N„ 



N -kN-k 2 N k 1+ k 2 -p N P 



N-kN-kA+k.-pNp 



Vk 2 N k 1+ k 2 -p N P VkA+k.-p^P 



U kl - U! k2 



'Jp-ki 



V-k 2 N k 1+ k 2 -p N P 



N k^-k 2 N kl+ k 2 -pNp 



}■ 



(2.10) 



Here, ui k = uj k (t), N k = N(k,t) and N k = 1 + N(k,t). 

It is instructive to notice that, by letting the finite initial time to — ► — co corresponding to the limit of complete 
collisions and using the identity sin(xt)/x ► ttS(x), Eq. (2.10) reduces to the standard on-shell Boltzmann 

t— >oo 

equation [19,30]: 



d 3 ki d 3 k 2 n5(u) kl + u> k2 -cj p _ fcl _ fc2 - u) p ) 

(2-7T) 3 (27r) 3 8u) kl UJ k2 UJp- kl - k2 

VkA+k^p 



N k ^k 2 N k 1+ k 2 -p N p 



}• 



(2.11) 



In this equation, the time integration in the original KB equation reduces to S function which represents the energy 
conservation in each binary collisions. Therefore, the quantum effect involved in the retardation is considered to be 
related to the energy broadening. It should be noted that while only the binary collisions contribute to the collision 
term in the on-shell Boltzmann equation because these collisions have no kinematical threshold on the mass-shell, 
<-> 4, 1 <-» 3 and 2^2 scattering processes are effective in the kinetic equations (2.9) and (2.10). 



D. Single particle excitation 

The kinetic equations obtained above are nonlinear equations for the Wigner distribution function. As an application 
to practical calculations, let us consider a single particle excitation: We add a particle with momentum p at t = 



4 



to a system which is in thermal equilibrium initially. We would like to calculate the relaxation rate for such an 
excitation. Since the distribution functions with momenta k p in the collision term don't change significantly 
from the thermal equilibrium value N(u> k ), we can take the self-energies (equivalently, the distribution functions with 
k / p) as in thermal equilibrium [30]. Under this situation, the kinetic equations reduce to linear equations for 
SN(p,t) (= N(p,t) — N(oj p j) as follows: Eq. (2.10) without the memory effect reduces to 

Q- t 5N{p, t) = -2 7 (p, t)6N(p, t) , (2.12) 

and Eq. (2.9) with the memory effect to 

g- t SN(p,t) = 2 J dz ^_ 7 (p,t- z )J 5N(p,z ) 

= -2j(p,t)SN(p,0)-2 [ dz 0l (p,t- z )^-SN(p,z ). (2.13) 

JO OZq 

Here, 7 (p, t) is a time-dependent damping rate given by 
A 2 f d 3 h d 3 k 2 1 



l{P, t) = 



u p J (2tt)3(2^ 



\2ojp J (2ir) 3 (2tt) 3 8uj kl uj k2 uj p - kl 



sinK+^+ Vfcl _ fa u p )t ^ M N M N {uJp _ ki _ k2) NMNMNi^^)] 



+ sinK +c fc +c p _ fcl _ fc +o, p )t w fe) -N(u, kl )N(uj k2 )N( o> p _ fel _ fe2 )] 

Wfej + Wfc 2 + cj p _ fel _ fe2 + w p 

+3 sinK +o; fc2 -o; p - fcl - fc2 -^ p )t [j^J^JJVK-fe-fa) _ AT (a;fe JAT(u, fe2 )^(u; p _ fel _ fe2 )] 
Wfej + ^fe 2 — W p _fei-fe 2 — ^ P 

^ sin^-i^-i^^-a^t ^ (wfeJ ^ KJ ^ (Wp _ fei _ feJ _^ KJ ^ (wfeJi v K _ fei _ fe2)]} ) (214) 
Wfe! - w fc2 - uj p - kl - k . 2 -uj p 

where distribution functions are in thermal equilibrium, i.e. N(u> k ) = l/(exp(u k /T) — 1) with temperature T, and 
uj k is independent of time. We have used the property j(p,t = 0) = from the first line to the second in (2.13). 
It it important that the right-hand side of (2.13) at some time t can be computed by integrating over the known 
functions for the time range z$ < t with given initial value of 5N(p, 0) thanks to 7 (p, 0) = 0. It is noted that this 
time-dependent damping rate is well defined even in a theory with long-range interactions like gauge theories in which 
self-energies in the on-shell limit are ill-defined [30]. 

In order to highlight the non-Markovian nature, it is instructive to consider the infinite-time limit t — > oo in which 
the time-dependent damping rate reduces to the on-shell damping rate 7oc (p) given by 

= _A^_ r rf 3 fci d 3 k 2 ir5(uj kl + uj k . 2 - uJ p - kl - k2 - u p ) 
loc[p} 4uj p J (2tt) 3 (2tt) 3 8uj kl u> k2 uj p - kl - k2 

x [N(cj kl )N(uj k2 )N(Ljp_ kl _ k2 ) - N(uj kl )N(Lj k2 )N(ujp_ kl _ k2 )] . (2.15) 

This quantity has been obtained analytically for p = [33] and numerically for p ^ [34,35]. In this limit, (2.12) 
and (2.13) coincide with the on-shell Boltzmann equation for a single particle excitation, 

^SN(p, t) = -2 loo (p)SN(p, t) , (2.16) 

whose solution represents an exponential damping of the excitation 

SN(p, t) = 5N(p, 0) cxp {-T-I(p) ■ t) , (2.17) 

with the relaxation time r re i(p) = l/2 7oo (p). The non-Markovian nature, or the memory effect, disappears in this 
limit. We, therefore, realize that it is related to the short-time phenomena. 
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III. NUMERICAL RESULTS 



In this section, we numerically calculate the thermal mass and the time-dependent damping rate for several values 
of the coupling constant. Then, we solve the kinetic equations with and without the memory effect. As expected, the 
relaxation time coincides that of the on-shell Boltzmann equation in the weak coupling regime due to a separation 
of scales. On the other hand, the results clearly demonstrate that the memory effect can influence relaxation of the 
single particle excitation in the moderately large and the strong coupling regimes. 



A. Renormalization of gap equation 

Since Green's functions in the self-energies are dressed ones as mentioned in Sec. II. A and the quasiparticle mass 
squared M 2 depends on the Green's functions, we have a a following gap equation for the thermal self-energies, 



M 2 = m 2 + Z a (x ,x ) 



2 A f d 3 k 2N(uj k ) + l 
= m + 



J (2tt)3 



2u>fc 



(3-1) 



where the thermal mass M is independent of time, and = \/k 2 + M 2 . 1 After eliminating the ultraviolet 
divergences in the momentum integral in the modified MS scheme, one obtains the renormalized gap equation which 
contains no UV divergences [36-38], 



M 2 = m 2 , 



\ r M 2 ( M 2 
In — 7T- 



1 



A,, f d 3 k N(uj k ) 



(27T) 



(3.2) 



Here, /i is an arbitrary subtraction scale, m r is the renormalized mass, and A r denotes the renormalized coupling 
constant which satisfies 



A~ 



din ^2 16-7T 2 



(3.3) 



This relation assures the solution M of (3.2) does not depend on /i. Below, we denote the renormalized mass and 
coupling as m and A for simplicity. By setting p, and m, M is obtained self-consistently for given A from Eq. (3.2). 



1 




g(27iT) 

FIG. 2. Thermal mass M over T is plotted as a function of the coupling constant. Renormalization scale is set to 2-kT. 



x For the thermal (static) self-energy, G < takes a same form in the KB and GKB ansatz, and consequently the corresponding 
gap equations are the same as (3.1). 



G 



Figure 2 shows the numerical results of M/T as a function of the coupling constant g with m = and p, — 2ttT. 
(Here and below, the coupling constant is rescaled via A = 24g 2 .) It is found that M/T increases monotonically as g 
is larger. This result agrees with the result in [36]. In the following calculations, we will use the solution M obtained 
here. 



B. Time-dependent damping rate 



Next we examine the time-dependent damping rate given by (2.14). Below, we focus on the zero mode p — 0. For 
p = 0, (2.14) reduces to 



7(P = 0,*) 
(24.g 2 ) 2 



1 



1 



dk\dk2 — 1 2 



| Si ((cj/^ + uj k2 + lo + - M) ■ t) - Si ((wfej + uj k2 + w_ - M) • t) 

-Si ((w fcl + w fe2 + u+ + M ) ■ t) + Si ((w fcl + cj fc2 + cj_ + M) • t) 

+3N{uj kl ) Si ((w fcl + uj k2 +oj+ + M) ■ t) - Si ((w fcl + w fe2 + w_ + M) • t) 

-Si {{uj kl + cj fe2 + cj+ + M) • t) + Si ((w fel + cjfc 2 + cj_ + M) • i) 
+Si ((— a;*! + Wfc 2 + lo+ - M) ■ t) - Si ((— Wfci + Wfc 2 + w_ - M) • t) 

+Si ((wfc 1 - ^fc 2 - w + — M) ■ t) — Si ((wfej - wfc 2 — w_ — M) • t) 

+3iV(w fcl )iV(wfe 2 ) Si((w fel +w fc2 + w+ +M) •<) - Si((cj fcl +^fe 2 +w- + M) -t) 

-Si((cj kl + uj k2 + M) ■ t) + Si ((w fcl + w fc2 + oj_ + M) ■ t) 

+2Si ((-ui kl + u k2 +uj + - M) ■ t) - 2Si ((-w fcl + cj fc2 + - M) • t) 
+2Si ((w fcl - io k2 M) ■ t) - 2Si ((uj kl - io k2 - W- - M) ■ t) 

+Si ((uj kl + uj k , 2 —u + — M) ■ t) - Si ((wfej + cj fc2 - w_ - M) • i) 

+Si((-w fcl - w fc2 + w+ - M) • t) - Si ((-wk, - w fc2 + w_ - M) ■ t) ] } , (3.4) 



where fci.2 = |fci,a|) UJ ± = \/ (&i ± fe) 2 + A/ 2 , and Si denotes the sine integral defined by Si(x) = Jq dtsxn(t)/t. 



g(27iT) = 0.1 




10' 



t*T 



FIG. 3. The time dependence of j(p = 0,i)/(24g 2 ) 2 T for three values of the coupling constant, g(2ivT) = 0.1 (solid), 
0.69 (dotted) and 1.61 (dashed). 

In Fig. 3, 7(p = 0,i)/24 2 <? 4 T is plotted as a function of t ■ T for three values of the coupling constant g(2,TrT) = 
0.1, 0.69 and 1.61. As we can see, for all values of the coupling constant, j(p = 0, t) oscillates in early time and then 
approaches to some constant value which corresponds to the value in the infinite-time limit, i.e. 7oo(p = 0). As the 
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coupling constant is larger, it is quicker for j(p — 0,t) to end oscillating. 2 



C. Decay of single particle excitation: comparison between the cases with and without memory effect 

Now let us examine the time evolution of the distribution function with/ without the memory effect. Below, we will 
call SN(p = Q,t) obtained from (2.12) "Markovian result' and that from (2.13) "non-Markovian result', and make a 
comparison between Markovian and non-Markovian results for several values of the coupling constant. 



1. Weak coupling regime 

Figure 4 shows the time evolutions of SN(p = 0,t) for the weak coupling g(2irT) = 0.1. As we can see, the 
Markovian and non-Markovian results agree completely. This is because a separation of scales is satisfied in the weak 
coupling regime: the relaxation time is much larger than the time at which the oscillation of the time-dependent 
damping rate ends. For g(2i:T) = 0.1, the relaxation time (times T) is about 4300, and oscillation of j(p = 0, t) ends 
at t ■ T ~ 20, as read directly from Figs. 4 and 3. In this case, SN(p — 0,i) essentially obeys the on-shell Boltzmann 
equation (2.16) and represents the exponential damping (2.17) in both cases. 




5000 ' ' ' mrooo 

t*T 

FIG. 4. 8N(p = 0, t)/SN(p — 0, 0) for Markovian and Non-Markovian results is plotted as a function oit-T. Non-Markovian 
result has been divided by 10. 



2. Moderately large coupling regime 

In Fig. 5, 8N(p = 0,t) /SN(p — 0,0) is plotted as a function of t ■ T in the moderately large coupling regime, 
g(2irT) — 0.69,0.91 and 1.02. In this regime, relaxation for the non-Markovian case is faster than that for the 
Markovian case, which is consistent with the results from the KB equations in 2+1 dimensions [19]. 



2 It should be noted that if we use the bare mass m instead of the thermal mass M in the sunset self-energies, the time 
evolution of j(p, t)/24 2 g 4 T is exactly the same for any value of the coupling constant. It is important to include a mean-field 
effect on the propagation of particles generated by the tadpole self-energy as well as a scattering effect inherent in the sunset 
self-energies on damping of the excitation. 
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t*T 

FIG. 5. Similar to Fig. 4 for the coupling constant g(2nT) = 0.69,0.91 and 1.02 from the top to the bottom. Solid and 
dotted lines represent Markovian and non-Markovian results. 



3. Strong coupling regime: Oscillation of the distribution function 

For the strong coupling constants g(2nT) = 1.04, 1.20 and 1.61, Fig. 6 shows absolute values of SN(p = 0, t)/SN(p = 
0, 0) as a function of t -T. We observe that SN(p = 0, t) for the non-Markovian case oscillates around zero for g > 1.04. 
3 As we can see in the middle panel, relaxation for Markovian and non-Markovian cases equals at g ~ 1.2. Then, 
relaxation for the non-Markovian case becomes slower than that for the Markovian case, as read from the bottom 
panel. 



3 We have checked that 5N with the memory effect doesn't oscillate for g < 1.04 by seeing the behaviors of 8N and dtSN 
in the late time. For g < 1.04, as time evolves, 8N/SN(t = 0) starting from unity and —d t 5N/SN(t = 0) from zero become 
parallel and never cross. They decrease monotonically in the late time, and 8N cannot reach the negative value. 
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FIG. 6. Absolute values of 8N(p — Q,t)/6N(p = 0,0) for Markovian (solid lines) and non-Markovian (dotted lines) are 
plotted as a function of t ■ T for the coupling constant g(2nT) — 1.04, 1.20 and 1.61, from the top to the bottom panels. 

D. Discussion 

Only in the weak coupling regime where damping of the excitation is exponential, the relaxation time r ro i can be 
obtained by an exponential fitting 

6N(p,t)~exp(-T^-t) . (3.5) 

In the moderately large and the strong coupling regimes, especially in the region where SN oscillates, we can not 
define the relaxation time of the zero mode in a precise sense. The important point, however, is that an opposing 
influence of the memory effect on relaxation of the excitation is observed as the coupling constant changes, as you 
can see in Figs. 5 and 6. 

Referring to the discussion in the weak coupling regime, it should be noted that the memory effect can affect the 
relaxation time even in the region where the separation of scales is satisfied well. For g = 0.69, although the relaxation 
time (~ 23 • T) in Fig. 5 is larger than ~ 3 • T at which the oscillation of j(p, t) ends in Fig. 3, there is the onset of 
the memory effect, that is, there is a slight difference between Markovian and non-Markovian results. Therefore, it is 
expected that the separation of scales works only in the adequately weak coupling. 

In Fig. 6, as the coupling constant increases, the relaxation time for the non-Markovian results seems inclined to 
saturate. On the other hand, the Markovian results show a monotone decreasing of the relaxation time. In order 
to see whether the relaxation time saturates, we need to study the time evolution for stronger coupling g > 1.6, 
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which is computable in the same way. However, such a strong coupling brings about the possibility of conflicting 
with the quasiparticle approximation used in our calculations, as mentioned below. Therefore, we don't make further 
consideration for the possibility of a saturation of the relaxation time here, which is an important future work. 

Throughout the present paper, we have used the quasiparticle approximation in which an excitation has an infinite 
lifetime. The validity of this approximation requires the condition that the relaxation time is larger than l/ui p 
(= l/M for zero mode), which corresponds to a long-lived excitation [30]. In the range of the coupling constant in 
our calculations, this condition arc satisfied. In order to examine the case of stronger coupling constant, it is needed 
to include the dynamical spectral function with the finite width, which is outside the present scope. 

IV. SUMMARY AND OUTLOOK 

Based on the kinetic equations with/without the memory effect, we studied the role of that effect in decay of the 
single particle excitation of 4> A theory. Starting from the two-particle irreducible effective action of the three-loop 
approximation and corresponding Kadanoff-Baym equations, we derived the kinetic equations into which the memory 
effect was incorporated by the GKB ansatz proposed by P. Lipavsky, V. Spicka, and B. Velicky [27]. It was observed 
in numerical results that, while relaxation of the excitation is unaffected by the memory effect in the weak coupling 
regime, this effect influences it differently as the coupling constant increases: The memory effect leads to faster 
relaxation for the moderately large coupling and to slower relaxation for g > 1.2. It is expected that the effect of 
memory manifests itself in totally noncquilibrium situation, that is, in non-linear kinetic equations for the distribution 
functions. The important point is that, although we adopted a simple (j) 4 model, the influence of the memory effect 
on relaxation is probably universal in strongly correlated plasmas. In particular, our results suggest a possibility of 
the saturation of the relaxation time. 

Finally, some comments are in order. The spectral function used in this paper is the quasiparticle one and cannot 
be applicable to the system where the relaxation time is shorter than 1/uj p . In order to study the effect of the finite 
width in the spectral function, it would be necessary to use the phenomcnological treatment such as an extended 
quasiparticle picture [32,39], or to solve directly KB equations for two-point Green's functions [13-16,18,19]. Two- 
time Green's functions in KB equations involve the information about the dynamical spectral function having the 
finite width as well as the one-particle distribution function. Furthermore, study of the linear sigma model or the 
fcrmionic system with the chiral phase transition is also an important future project in our approach, which might 
be related to Heavy-ion collision experiments. 
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